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Abstract 

We continue to study a cosmological model with large-scale inhomo¬ 
geneity. Working in the cosmic rest frame we determine null geodesics, 
redshift and area and luminosity distances. Combining the result with 
Hubble’s law enables us to calculate the distance of the local group of 
galaxies from the origin r = 0 where the Big Bang has taken place. 
We obtain a surprisingly small value of about 2 million lightyears. 
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1 Introduction 


In this paper we continue the study of inhomogeneous cosmology in the 
cosmic rest frame. The perturbative solution of Einstein’s equations in¬ 
vestigated in [1] is now crucially tested by determining the corresponding 
redshift - distance relation. We work in the cosmic rest frame throughout. 
The big advantage of this frame over comoving coordinates is that it is really 
a frame, that means a global coordinate system for the whole universe. As 
a consequence its origin r = 0 is the point where the Big Bang has taken 
place, and t = 0 is the time 14 billion years ago when this has happened. A 
natural and meaningful question which now occurs is: Where is r = 0 with 
respect to our place in the Milky Way ? In [1] this question remained open. 
Here we will answer it: Surprisingly enough r = 0 is not further away than 
the Andromeda galaxy ! This follows by taken Hubble’s law into account in 
the inhomogeneous cosmological model. 

The paper is organized as follows. In the next section we recall the main 
results of [1] and use it to determine null geodesics; these are needed to 
find the light rays in the expanding cosmic gravitational field. In section 3 
we consider Maxwell’s equations in the eikonal approximation and calculate 
the redshift. In section 4 we compute the area and luminosity distances for 
an observer near the origin r = 0. Combining the result with Hubble’s law 
enables us to determine the distance R of the observer from the origin. We 
obtain a surprisingly small value for R of about 2 million lightyears. 

2 Null geodesics 

We choose spherical coordinates t, r, •&, (f> in the cosmic rest frame and assume 
the inhomogeneous metric of the form 

ds 2 = dt 2 + 2 b(t, r)dtdr — a 2 (t, r)[dr 2 + r 2 (<A9 2 + sin 2 fide/) 2 )]. (2.1) 

In [3] we have constructed a perturbative solution of Einstein’s equations 
given by 

a(t, r ) = at 2 / 3 + tg\(r) + h\(r) + 0(a -1 ) (2.2) 

b(t,r) = at 2 / 3 Lr + tg 2 (r) + h, 2 {r ) + 0{oT l ). (2.3) 

Here a and L are universal constants of integration and the functions g \, g 2 
and hi, h 2 are explicitly known. For later use we also list the components 
of the metric tensor corresponding to the line element (2.1) 

3oo = 1, 3oi = b(t,r), g n = ~a 2 (t,r) 
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(2.4) 


922 = ~r 2 a 2 (t, r ), g 33 = -r 2 a 2 (t , r) sin 2 d, 

and zero otherwise. The components of the inverse metric are equal to 


„„ _ n _ 1 

9 ~ D , 9 9 ~ D 


„00 


.01 


9 22 = 


a 33 = 

2r.2 ’ y 


1 


err 


a 2 r 2 sin 2 x) 


where 


D = a 2 + b 2 


(2.5) 

( 2 . 6 ) 

(2.7) 


is the determinant of the 2x2 matrix of the t, r components. The non¬ 
vanishing Christoffel symbols are given by 


r{jo = -|, r° 1 = -^a, r° n = ±(a 3 a-aba' + a 2 b') (2.8) 


T ^2 = — (ra 3 a + ba 2 + rbaa'), T^ = ———(ra 3 a + ba 2 + rbaa' ') 

r oo = r oi = 5 rjj = -^(a&d + aa' + W/) 

T ^2 = — (raba — a 2 — raa'), T^ = — (raba — a 2 — raa') sin 2 d, (2.9) 


Tq 2 = T 2 2 = — + r §3 = — sin d cos d 

u a a r 33 

r 3 _ a r 3 _ a ' , 1 r 3 _ cos ^ 
io 3 “a’ il 3 “a + r’ 123 “ ii^T 


( 2 . 10 ) 

( 2 . 11 ) 


Here the dot means d/dt and the prime d/dr. 

The propagation of radiation in the cosmic gravitational field is described 
by the geodesic equation 


7 + 


where 


AT = 


dx y 
dr 


( 2 . 12 ) 


(2.13) 


is the wave vector and r is the affine parameter along the ray. We want to 
determine the null geodesics which starts from a point on the z-axis under 
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a certain angle and lies in a plane (f) =const. Then only k°,k 2 and k 2 are 
different from zero and the geodesic equation in leading order read as follows 


dk° L 2 r 2 d,, n , 9 Lr , n oa.,,., 

~f + ~ a {k> 


where 


+—(rad + Lr)(k 2 ) 2 = 0 
P 

dk Lv d , r \ x 9 ^ a T n , i 1 _ . 9 \ / , 1 x o 

~F - F^ {k ) + V + t ra + L r)(k ) + 


(2.15) 


+ 1 —(Lr 2 ad — 1 ){k 2 ) 2 = 0 

P 

(2.16) 

dk a j a o 2 i o 

— + 2—k°k 2 + -k 2 k 2 = 0 
or a r 

(2.17) 

P = L 2 r 2 + 1. 

(2.18) 

immediately be integrated 


2 C dtf 

(ar ) 2 dr 

(2.19) 


where C is a constant of integration which fixes the angle under which the 
trajectory starts at the z-axis. In the following we are interested in small 
C, that means the trajectory is nearly radial. To solve the two remaining 
equations (2.15-16) we use the ansatz 


k° = — + fo(t,r), k 1 = K — + fi{t,r). (2.20) 

a a- 

Here K is another integration constant which multiplies the radial null 
geodesics. The latter follows easily from 

ds 2 = 0 , kijk/ 1 = 0 , ( 2 . 21 ) 


it corresponds to C = 0. The sign of the square root has been chosen 
in such a way that the light ray moves towards r = 0. For small C we 
neglect quadratic term in fj, then we obtain from (2.15) the following linear 
equation 


dfo of KLr d 0 K d _ C 2 d 
dr 10 y//3 a 2 Jl y/f3a a 3 /3r 2 


( 2 . 22 ) 
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The normalization (2.21) gives another linear equation between fo and f\: 


2Kf iy /p = -^-2f 0 -{P-Lry/p). (2.23) 

a z r z a 

Substituting this into (2.22) we find the simple equation 

^ + 2/ 0 ^ = 0 (2.24) 

cLt a z 

with the solution 

fo = (2-25) 

a z 

where I\\ is another constant of integration. 

Now by (2.23) we get the second function f\ 

fl = 2Ka 2 r 2 y/~j3 (2 ' 26) 

Since 

k o = dt = h +f ° (2.27) 

dr a 

the constant K\ in (2.25) only enters in the dependence f(r). Therefore, 
without loss of generality we can put K\ = 0 by a redefinition of the affine 
parameter. Finally we have to check that the second geodesic equation (2.16) 
is satisfied by (2.25-26). This is indeed true in leading order 0(a -3 ), higher 
orders do not interest us because we have neglected them in the Christoffel 
symbols already. 

Summing up the wave vector has the following components 


We put 


k° 


I< 

a 


dt 

dr 


1 £r - va <T dr 

a 2 2 Ka 2 r 2 -\[/3 dr 

2 _ _C_ _ cM 
a 2 r 2 dr' 


C 


ro= K 


(2.28) 

(2.29) 


(2.30) 
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and assume this to be positive; its physical meaning will soon become clear. 
By dividing (2.31) and (2.28) we get 


dt) ro 

dt ar 2 


(2.31) 


Since this is positive the polar angle 'd(t) is monotonously increasing with 
time. We consider a trajectory starting on the z-axis at a point r\ at time 
t\ so that d{t\) = 0. This is the place of our radiating galaxy. To find the 
light path we divide (2.29) by (2.28) 


¥ 

¥ 




dt 


2 r 2 yf]3J 


(2.32) 


and separate the variables: 


dt 


r 2 VL 2 r 2 + 1 


-dr. 


(2.33) 


a(t) Lr^yJ L 2 r 2 + 1 — r 2 (L 2 r 2 + 1) + Tq/2 

. We consider the light rays towards the origin r = 0, that means 

dr 

— < 0 . 

dt 

Setting this equal to 0 we find from (2.32) the minimal distance of the path 
from the origin 


4 2 

r 0 , r 0 


= y + Vt + p- 


(2.34) 


so ?’o determines this minimal distance. 

To find the path we must integrate (2.33). To make this integration 
simple we assume 

(2.35) 


lV > l 


and expand 


\ffd = Lr + 


1 


2 Lr 


The usefulness of this approximation was discussed in [1], We integrate 
(2.33) from an emission time t,\ to the present time T: 


= —(T 1/3 -t{ /3 ) = -2 L 


dt 


at?/ 3 3a 


R 


1/3 s 


r* — 7n 


>dr 
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(2.36) 


Here r\ and R are the radial coordinates of the emitter and observer, re¬ 
spectively. Solving for t\ we obtain the equation t = t(r ) of the light path 

t = \T 1/3 - 3 aL{r 2 - R 2 + rg log Z~ T \ )1 3 . (2.37) 

3 The redshift 

The propagation of electromagnetic radiation in the expanding universe is 
described by Maxwell’s equations in a gravitational field. The homogeneous 
equations are 

V \F /lL/ + V fjF v x + V v F\ t j = 0 (3.1) 

where V stands for the covariant derivative with respect to our inhomoge¬ 
neous metric. In addition the field tensor F^ u satisfies the inhomogeneous 
equations 

= 0 (3.2) 

away from the radiating galaxy. As in Minkowski vacuum, equation (3.1) 
implies the existence of a vector potential 

Ffu, = V m A, - (3.3) 

and similarly for the upper indices. Here we can impose the “Lorentz” gauge 
condition 

= 0. (3.4) 

As in classical electrodynamics we write (3.2) in terms of the vector 
potential 

= 0. (3.5) 

Now we want to commute the derivatives in the second term and use (3.4). 
Since the covariant derivatives do not commute we obtain an additional term 
with the Ricci tensor 

- R"A a = 0. 

In cosmology it is sufficient to solve this equation in the eikonal approxima¬ 
tion by considering a solution of the form 

A u (x) = f u (x)sin[(p(x)}. (3.6) 
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Here f v is a slowly varying amplitude, p is the phase which varies from 0 
tO 2ir over the wave length of the radiation. Since the latter is extremely 
small compared to cosmic distances, tp is big so that we shall consider only 
the first two orders of ip in (3.6). The leading order 0(p 2 ) gives 

V^VV^ = 0 (3.7) 


and in order 0(<p ) we obtain 


2 V^rv^ip + f'V^V'V = 0. (3.8) 

As in optics one introduces the wave vector 

= VV (3.9) 

Then the two equations to be solved read 


o 

II 

3. 

(3.10) 


(3.11) 

From equation (3.10) we have 


WVvkp = 0, 

(3.12) 


and since (3.9) implies 

V lykfi = V 


we arrive at 


k! t V lt k v = 


(3.13) 

Let now x M (r) be the curve 


dx>1 - u> 

~d^~ k 1 

(3.14) 

then we have 

dk v 

d 2 x u dr d 2 x v 1 


dx^ 

dr 2 dxv dr 2 kJ‘ 


and (3.13) becomes 

d 2 x u 

dx^ dx x 

(3.15) 

dr 2 

+ » x dr dr ~ °' 

This is the geodesic equation (2.12) for x v (t) which we 

have solved in the 


last section. 



Following G.F.R. Ellis [2] (see also [3]) the redshift is given by 


^em _ i _ (^^)em 

Robs WV obs 


(3.16) 


where 

tively. 

(6.9): 


is the 4-velocity of the emitting galaxy and the observer, respec- 
The radial velocity of the galaxy we know from [1] equations (6.8)- 


uq = sf]3, u\ 


L 2 rt 

w 


(3.17) 


Then, since 


u = u^kP = K —7^r [ 1 + 
a(t) L 


L 2 rt Lr — VP' 

P a(t) 


(3.18) 


the redshift becomes 


1 + z = 


VP O(T) 


VPr a(t) L p a (t) 


1 + 


L 2 rt / 


(Lr - VP + 


2r 2 VP 


x 


1 + 


l 2 rt 

Prcl(T) 


(/LR — sfpR + 



2 R 2 Vfa 



(3.19) 


where T is the present age of the Universe and R > 0 the radial coordinate 
of the observer. The emission time t follows from the coordinate distance 
r = r\ of the galaxy by means of (2.37). For later use we introduce the 
expression 

*, = 1 + ^ (ir _^ + _rL_). (wo. 

Then the redshift is equal to 



, _ VP a(T) K r 

VPr a V) Kr 

(3.21) 

where 

Pr = L 2 R 2 + 1 

(3.22) 

according to (2.18). 
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4 Area and luminosity distances 


In this section we must specify the radial coordinate R of the observer. We 
choose 

R = r ' m ; n (4-1) 

equal to the minimal distance (2.34) of the light path from the origin r = 0 
where the Big Bang has taken place. This greatly simplifies the following 
calculations because we have 


dr 

dt 

= 0 

R 

(4.2) 

Kr = 1 

(4.3) 

so that the redshift (3.21) becomes 



i +2 = 

VPr a(t ) 

(4.4) 


Theoretically the simplest radial distance is the area or angular-diameter 
distance [4] 

D a = a(t)r(t ) (4.5) 

where r(t) is the radial coordinate of the galaxy at the time of emission 
t which can be found by inverting (2.37). Assuming reciprocity [3] Da is 
related to the luminosity distance Dl which is measured by the astronomers 
by 

D l = (1 + z) 2 D a = a(t)r(t)( 1 + z) 2 . (4.6) 

Considering t as a parameter the equations (4.4) and (4.6) determine the 
distance redshift relation Dl = Dl{z ) in parametric form. 

We want to calculate the derivative 


(ID /\ _ c 

dz 2=0 Hq 


(4.7) 


which gives the Hubble constant Hq at present time T (we have put the 
velocity of light c equal to 1). From (4.6) we obtain 


dD L 


dz 


= ci (T) C 4- R + 2a(T)R 

2=0 v dz T K ’ 


(4.8) 


where (4.2) has been taken into account. Again using (4.2) we get from (4.4) 


dz dz 
dt T dt 


a 

a 


+d t K r | T . 

T 


(4.9) 
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Since 


dtK r \r = 0 

because r(T) = R is the minimal distance (4.1) we get in leading order 0(a ) 


( 2 . 2 ) 

dz 

2 

(4.10) 


dt 

T 3 T' 

Using 


o „ 



■ <rr\ z u 
«( ) 3T 

(4.11) 

we finally find 

dD L 

dz 

= a(T)R. 

z=0 

(4.12) 

Now it is time to put some simple numbers in. Let 



£?crit = 1.878 x 10 29 h 2 g/cm 3 

(4.13) 

be the critical density [4] and we assume a Hubble constant 



h = 0.7, Hq = 70 kms : Mpc 1 . 

(4.14) 


We use the empirical fact that the corresponding Hubble time 

—— = 14 x 10 9 years = T (4.15) 

Hq 

coincides with the age T of the universe. A realistic density of normal matter 
at present time T is 


g m - 0.01 x £ crit - 67tG/3rT 2 ( 4 - 16 ) 

where the last equality is the lowest order result [1] eq.(6.7). This allows to 
determine 


Pr = L 2 R 2 + 1 = 44.3 (4.17) 

which gives 


LR = 6 . 6 . (4.18) 

Next we assume a radial velocity of our Galaxy of 300 km/sec that means 

LR 

Vm = 0.001 = —- 4.19) 

a{T) 
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where the last equality comes from [1] eq.(7.6). Using (4.18) this allows to 
determine 

a(T) = 6.6 x 10 3 . (4.20) 

Now we are able to calculate R from (4.12) using (4.7) (4.12) and the 
Hubble time (4.15) 


T 

R = —— = 2.1 x 10 6 ly = 640 kpc. (4.21) 

This is the distance between our Galaxy and the origin r = 0 where the Big 
Bang has taken place. It is even smaller than the distance of the Andromeda 
galaxy which is 2.5 x 10 6 lightyears away. From the dipole anisotropy of the 
CMB one has derived a net velocity of the local group of galaxies of 627 ± 22 
km/sec in a direction between the Hydra and Centaurus clusters of galaxies 

[4] [6]. We expect that the origin r = 0 lies in the opposite direction. The 
velocity about 600 km/sec seems to be rather high. But this velocity is the 
sum of the cosmic radial velocity and a local velocity due to the attraction of 
the Hydra and Centaurus clusters. Since the latter has not been estimated 
we have made the simple assumption that 300 km/sec is the cosmic velocity 
(4.19). 

If we calculate the second derivative c^D^/dz 2 in the same way we get 
zero. This shows that a lowest order calculation is not good enough to 
determine the acceleration parameter t/o- 
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